Self-Diffusiophoresis in the Advection Dominated Regime 
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In both biological and artificial systems, concentration gradients can serve as a convenient mech- 
anism for manipulating particles and generating motility. Particles that interact with a solute will 
move along its gradient; if they themselves generate the gradient, this mechanism provides a means 
of self-propulsion. We consider a version of this type of motility appropriate to certain biological 
systems where polymeric filaments provide the concentration gradient. As the filament diffusion is 
small, this corresponds to a regime of large Peclet number where the motion is dominated by the 
effects of fluid advection. The nature of such concentration-gradient-driven motion in the advective 
regime differs in certain fundamental respects from the same process at low Peclet number. In 
particular, we show that out of four broad scenarios of steady state motion at low Peclet number, 
only two remain viable in the strongly advecting limit. 

PACS numbers: 87.19.ru, 47.63.-b, 82.70.Dd, 47.10.-g 



Certain motile biological objects, such as the bac- 
terium Listeria monocytogenes and replicating chromo- 
somes in asymmetric bacteria such as Caulobacter cres- 
centus and Vibrio cholerae, generate their propulsion by 
polymerizing or depolymerizing protein filaments [iHS]. 
This establishes a concentration gradient in the protein 
so that there are more filaments on one side of the ob- 
ject than the other. In the case of Listeria, the activa- 
tion of the Arp2/3 protein complex at the cell posterior 
promotes the polymerization of actin and results in the 
formation of a distinctive comet tail, a region of high 
actin concentration immediately behind the moving bac- 
terium [T]. Conversely, in Caulobacter and Vibrio the 
chromosome stimulates the disassembly of filaments of 
ParA so that the concentration of ParA is higher in front 
than it is behind the chromosome as it translocates |4j. 

A colloid that maintains an asymmetric concentration 
of solute around it will propel itself through a fluid with 
a well-defined velocity. This phenomenon, known as self- 
diffusiophoresis, has been exploited as a means of propul- 
sion of micro- or nano-swimmers [5"l2|. In a typical 
example, the concentration gradient is controlled by an 
active region on the colloid that catalyzes a chemical re- 
action, leading to more product and less reactant near 
the active region than on the far side of the colloid. An 
interaction between the colloid and solute sets up fluid 
flow within a thin boundary layer with a slip velocity 
at the edge of this layer proportional to the concentra- 
tion gradient [71 |T3] . This local concentration gradient 
is in turn determined by a combination of diffusive and 
advective fluxes, whose relative importance is measured 
by the Peclet number Pe = Ua/D, where U and a are 
a characteristic velocity and length scale of the colloid, 
respectively, and D is the solute diffusion constant. In 
experimental realizations considered to date (such as the 
decomposition of H2O2 [8 J, the solute has a large diffu- 
sion constant and Pe is small (^ 10^"^) so that the con- 
centration gradient is effectively set by diffusion alone. 



Similarly, theoretical treatments have considered only the 
limit Pe=0 [71 [SHIl] • In this limit it is known that only 
two conditions are needed to achieve motility: the motile 
object must be able to maintain an asymmetric solute 
distribution in steady state, and there must be a net in- 
teraction between the solute and the object. However, in 
the biological examples of interest, the relevant solutes 
are filamentous protein assemblies with effectively van- 
ishing diffusion constants and it is not known whether 
the same conditions apply for motility. 

In this Letter, we determine the general properties of 
propulsion in a self-generated concentration gradient in 
the limit of infinite Peclet number. We find that ad- 
vective transport through the boundary layer leads to 
significant differences in the conditions needed to sustain 
steady state motion and in the scaling of the speed with 
surface activity, as compared to the Pe=0 limit. Some 
of our results are consistent with known observations of 
polymerization- or depolymerization-driven motility in- 
volving filamentous proteins, and the remainder can po- 
tentially be tested in such systems. 

We consider the idealized minimal model of Ref. [7], 
in which a spherical particle of radius a produces solute 
from an active patch on its surface at a rate a, shown 
schematically in Fig. [T] We restrict our attention to ax- 
isymmetric, steady state motion, viewed from the rest 
frame of the particle. The particle interacts with the so- 
lute through a potential V{r) with a range 5 -^ a that 
is small compared to the particle size; this separation of 
length scales allows the effect of the interaction to be de- 
termined through a boundary layer analysis that neglects 
surface curvature on the scale S. Consequently, the flow 
within the boundary layer is well approximated by the 
solution of the Stokes equations 



= -Vp + ^V^u - cVF, V • u = 0, 



(1) 



where fj, is the fluid viscosity and c the solute concentra- 
tion, in the upper half space (z > 0) that approximates 
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FIG. 1: (Color online) Schematic illustration of self- 
difFusiophoresis for a colloid at low (left) and high (right) 
Peclet number indicating the form of the slip velocity. The 
magnified region shows the nature of the approximation used 
in the boundary layer analysis. 



the local environment near the particle surface, as seen 
on the scale 5. The resultant flow is tangential to the 
surface with a limiting value, known as the slip velocity, 
at the outer edge of the boundary layer of 
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dz-z^V||c(z,x)a,V^(z), (2) 



where Vn denotes gradients projected into the tangent 
plane of the surface. When Pe is small and the solute 
diffuses rapidly it is a good approximation to assume 
that the concentration maintains a local equilibrium with 
a Boltzmann distribution c{z,x) = c(<5, a;) e"^*^^)/*^^-^, 
which reduces the slip velocity to u^^'^p = TO^Vyc where 
TO^ = ^i£^ j^Azz[l - Q-vi^)/kBT^^ is a diffusiophoretic 
mobility [TJ]. 

At arbitrary Pe, the solute concentration field obeys a 
conservation equation outside the boundary layer 



<9tc + u • Vc - I^V^c = 0, 
together with the boundary condition 



luf^c - 



DdrC\ = 



r— a+(5 



(3) 



(4) 



expressing that in steady state the rate of production 
of solute at the particle surface must equal the rate at 
which it is transported across the boundary layer and into 
the bulk fluid. At small Pe, the advective terms in Eqs. 
(Is]) and ^ may be neglected, while for thin boundary 
layers the integral boundary condition is replaced with a 
pointwise equality —DdrC = a j^, which yields an unique 
solution for the concentration field. In the opposite limit 
of large Pe the solute is still conserved and Eqs. ([s]) 
and (l4| still apply. However, if diffusion is negligible the 



solute can only be transported away from the particle 
surface if it is advected by a radial flow. This radial slip 
velocity, neglected in the tangential projection of Eq. ([2| , 
is therefore crucial to retain in the advective regime. 

However, the presence of a radial slip does not invali- 
date Eq. ([2]), as this is only an asymptotic formula that 
neglects surface curvature on the scale S. Dimensional 
analysis of the fluid continuity equation shows that the 
radial slip is typically smaller than the tangential slip by 
a factor S/a, so that its neglect is based on the short range 
of the interaction of the solute with the particle surface. 
Since this range remains short at high Pe, Eq. (pi) will 
remain a good approximation, with one caveat, namely, 
that any tangential flow around a spherical body must 
vanish somewhere (and the sum of the indices of all the 
zeros must be 2 - Poincarc-Hopf theorem) and in the 
vicinity of any such point the neglect of the radial slip is 
not justified. 

Away from these zeros the slip will continue to be ap- 
proximated by the tangential boundary layer flow, which 
should behave as in Eq. [2] 



slip 
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a + S 



dec{a + 6,9), 



(5) 



in spherical coordinates, where m is an 'advective mo- 
bility.' A simple but significant observation coming from 
this relation is that there is no tangential flow where there 
are no gradients. At low Pe diffusion spreads the solute 
throughout space, at least close to the particle. However, 
if the solute does not diffuse then it will only be found 
where it is produced, or has been transported by the fluid. 
Thus if we consider an active patch, covering the portion 
< 6 < 6p oi the sphere, the vanishing of the solute con- 
centration everywhere outside the patch implies that the 
tangential slip velocity also vanishes there. Fluid conti- 
nuity then requires that the tangential flow within the 
patch region is supplied by a radial influx at the edge 
9 — 6p and subsequently escapes via a radial outflux lo- 
calized near some internal point, 9 — hy axisymmetry. 
These considerations determine the qualitative features 
of the flow at high Peclet number. 

To be more quantitative we take this picture as an ide- 
alization, so that there is radial flow only at 9 = 9p and 
9 = 0, and purely tangential slip elsewhere. Fluid con- 
tinuity constrains this tangential flow so that sin(0)wg '^ 
is a constant, which we write as —A with A > 0. Eq. 
(Is]) then allows us to solve for the solute concentration 
profile at the edge of the boundary layer 



c{a + S, 



X{a + S) 



In 



tan{^9p) 
tan{^9) 



(6) 



The singular behavior as — >■ is an artifact of neglect- 
ing the radial slip. We regularize this by cutting off the 
solution at some small angle 9q within which the concen- 
tration assumes a constant value. 



The net radial outflux from the boundary layer is equal 
to the net tangential influx across a cylindrical surface 
spanning the thickness of the boundary layer at any 
6q < 9 < 9p, which for thin boundary layers we approxi- 
mate as —2'KaSPsm{9)uy^, where /3 can be expected to 
be an 0(1) constant coming from the average over the 
boundary layer. Thus the radial fluid outflux is given by 
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(7) 



and is small compared to the tangential flow by a factor 
5 /a as expected. 

At low Pe, the radial flow from the boundary layer can 
be neglected because 5 /a <^ 1. However, it is a feature of 
the high Pe regime that this radial outflux provides the 
sole mechanism for the transport of solute produced at 
the surface through the boundary layer and into the bulk 
fluid. This transport is given by Eq. ^ with diffusion 
neglected and leads, for uniform activity a{9) = a within 
the active patch < < 6'p, to an expression for the 
tangential flux A 



A^ = 



2am^asm^{^9p) 

(3{a + S)S 
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Finally, the speed of the colloid - moving along the neg- 
ative z-direction - as well as the entire flow field outside 
the boundary layer may be obtained from Lighthill and 
Blake's analyses of squirming spheres P^[T3] : 



U 



Asin2(i0p), 



(9) 



where we have neglected contributions from the radial 
slip that are smaller by a factor of 0{6/a), and the fluid 
flow field outside the boundary layer is shown for a Janus 
particle (6'p = f ) in Fig. [5] 

Tacit in the foregoing analysis were the assumptions 
that the particle was acting as a source for the solute, 
a > 0, and that the interaction potential was repulsive, 
V^ > 0. Reversing these signs independently provides four 
broad classes of self-diffusiophoretic motion, all of which 
lead to steady state motion at Pc = [7]. However, at 
Pe = oo only two cases continue to support steady state 
propulsion; repulsive producers (F > 0, a > 0) and at- 
tractive consumers {V < 0, a < 0). The reason for this 
again lies in the advective transport through the bound- 
ary layer in which radial slip transports material away 
from, or towards, the surface according to whether the 
activity is positive (producing) or negative (consuming). 
By fluid continuity any such radial flux will be converted 
into a tangential one with m^ '^ negative if there is an out- 
flux and positive if there is an influx. Thus the particle 
motion will always be away from solute that it is produc- 
ing, and towards solute that it is consuming, regardless of 
the interaction. But this motion is in the same direction 
as predicted by the boundary layer analysis of Eq. (pi) 




FIG. 2: (Color online) The flow field of a self-diffusiophoretic 
Janus particle at small (top) and large (bottom) Pe. The 
flow is shown scaled in units of the translational speed U and 
distances are scaled in units of the particle radius a. Note 
that the flow exhibits a fore-aft symmetry at low Pe, which is 
broken when the Peclet number becomes large. 



only if the interaction is repulsive in the former case and 
attractive in the latter, so that only these cases provide 
consistent steady states at high Pe. 

The conditions for motility in the Pe = cxd limit are 
therefore more restrictive than those at Pe = 0. It is no- 
table that the biological phenomena mentioned earlier, 
namely actin-driven propulsion of Listeria and chromo- 
somal translocation in Caulobacter and Vibrio, appear 
to coincide with the two cases that lead to steady-state 
motility in the high Pe limit. In Caulobacter, for ex- 
ample, the chromosome attaches to ParA via parS and 
ParB |4j, suggesting that the interaction between the 
chromosome and ParA is attractive, so that the chro- 
mosome acts as an attractive consumer of ParA fila- 
ments [T5] . In the case of Listeria, the bacterium surface 
and actin are negatively charged, so there should be a 
screened Coulomb repulsion between the two. While the 
actin comet tail is attached to the bacterium, it is plausi- 
ble that binding is sufficiently rare and transient [T71 [TS] 
that there is a well-defined net interaction that is re- 
pulsive, making Listeria a repulsive producer of actin 
filaments. We note that other mechanisms have been 
advanced for the case of Listeria motility |19j in addi- 
tion to self-diffusiophoresis [20] • Since even at high Pe, 
asymmetric solute distributions give rise to motility with 
appropriate signs of interactions, the mechanism of self- 



diffusiophoresis must be reckoned with. 

The speed ^ found at high Pe should be compared 
with the corresponding expression 



U 



AD 



sin2(0p) 



(10) 



for an identical active patch and interaction at low Pe [7] . 
Two observations are particularly noteworthy. 

First, the two limits exhibit distinct dependences on 
the size of the active patch Op. In the diffusive regime 
the speed is invariant under Op ^ t: — Op reflecting the 
fact that it is really the asymmetry that controls motility 
and, in particular, that the particle does not move if the 
coverage is uniform. At high Pe the situation is different, 
although some care is needed to interpret Eq. (Is]) for the 
tangential flux constant A. It is appropriate to rewrite 
the logarithmic term as ln(a/(5), which exhibits the rel- 
evant scaling without the artificial singularity associated 
with the neglect of the radial slip near 6 = tt. With this 
moUifier the speed in the advective regime has the depen- 
dence U ^ sin^(^6'p) and thus does not vanish in the limit 
0p — > TT of complete coverage. One interpretation for this 
result is that a uniformly active particle, which has no 
natural directionality, is unstable in the limit Pe = oo 
with respect to a spontaneous symmetry breaking that 
selects a direction for motion. Such spontaneous symme- 
try breaking has been observed in the case of actin-driven 
motility if the colloid radius is small enough [H [21] . This 
qualitative agreement with the present analysis encour- 
ages a more thorough analysis of the onset of motility to 
explain the observed size dependence. 

Second, the two expressions exhibit different scalings 
with the activity rate a; in the diffusive regime, U ^ a, 
while in the advective regime, U ~ a^". This discrep- 
ancy reflects precisely the different fluxes in Eq. Q. 
The diffusive flux is proportional to the tangential slip, 
whereas the advective flux is proportional to its square, 
and since the speed is proportional to the tangential flow 
the scalings follow. Recent experiments show that the 
propulsion speed of Listeria depends on the distribution 
of Act A, which catalyzes actin polymerization [52]; the 
speed increases with the amount of ActA but it is diffi- 
cult to adduce a functional form for the dependence of U 
on a from the data. 

In summary, we have shown that objects that are 
repulsive producers or attractive consumers of non- 
diffusing solutes should propel themselves through a 
fluid. Our abstract, minimal analysis suggests that in 
such situations, objects can spontaneously break symme- 
try to propel themselves in a given direction and achieve 
sustained steady-state motility. Note that we have not 
taken into account the effect of the non-diffusing solute 
on the flow; for example, the actin comet tail of Liste- 
ria should behave as a porous medium to suppress flow 
within it. Since this effect could be modeled with an in- 



ert trailing particle, which does not prevent motility at 
low Pe [22, it seems unlikely to be a significant factor 
at high Pe, but this avenue should be pursued in future 
work. Another important issue is the stall force needed 
to prevent motion. At low Pe, the stall force is simply the 
drag force, but at high Pe it will be determined by the 
distortion of the concentration profile by the boundary 
layer flow. 
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